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Abstract 

We study the entanglement entropy in confining theories with gravity duals using the 
holographic prescription of Ryu and Takayanagi. The entanglement entropy between a 
region and its complement is proportional to the minimal area of a bulk hypersurface 
ending on their border. We consider a disk in 2+1 dimensions and a ball in 3+1 dimensions 
and find in both cases two types of bulk hypersurfaces with different topology, similar to 
the case of the slab geometry considered by Klebanov, Kutasov and Murugan. Depending 
on the value of the radius, one or the other type of hypersurfaces dominates the calculation 
of entanglement entropy In 2+1 dimensions a useful measure of topological order of the 
ground state is the topological entanglement entropy, which is defined to be the constant 
term in the entanglement entropy of a disk in the limit of large radius. We compute 
this quantity and find that it vanishes for confining gauge theory, in accord with our 
expectations. In 3+1 dimensions the analogous quantity is shown to be generically nonzero 
and cutoff-dependent. 
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1. Introduction and summary 

Certain systems in two spacial dimensions are known to exhibit topological order [1,2]. 
In this case the ground state of the system has highly non-trivial properties, related to the 
degeneracy on higher genus surfaces and the existence of quasi-particle excitations. Such 
a non-local order can exist even when local order parameters vanish. Recently a quantity 
called "topological entanglement entropy" was proposed [3,4] to measure the degree of 
topological entanglement. It involves computing the entanglement entropy of a large disk 
and extracting the constant term, which does not scale with the circumference of the disk. 
The natural generalization to higher dimensions would involve entanglement entropy of a 
large ball. 

Computing the entanglement entropy for a general system is not a trivial exercise. 
Recently Ryu and Takayanagi [5] proposed a way to compute the entanglement entropy 
in theories with holographic dual. The computation involves finding a certain minimal 
surface in the bulk of the asymptotically AdS space. 1 

This hypersurface approaches the border between the region and its complement on 
the boundary. The generalization of this prescription also exists for the theories whose 
duals are not asymptotically AdS. In this paper we make use of this proposal to compute 
the topological entanglement entropy of some simple confining theories in 2+1 and 3+1 
dimensions. 

More precisely, we will study D3 and D4 branes compactified on a circle, which at 
large distances become confining 2+1 and 3+1 dimensional theories respectively. The 
entanglement entropy of the "slab" (the subspace defined by — 1/2 < x < £/2, where 
x is one of the spacial coordinates) in these systems has been studied in [15]. As it 
turns out, there are two types of minimal surfaces: connected and disconnected ones. In 
ref. [15] a first order phase transition between these two types of solutions was observed 
as a function of £. For small I the connected solution dominates the computation of 
entanglement entropy, while for large I the disconnected solution becomes preferred. The 

1 Schwimmer and Theisen [6] argued that anomalous terms in the entanglement entropy are 
not correctly reproduced in the holographic prescription. However the validity of this claim has 
been questioned in [7]. The strong subadditivity of the holographic entanglement entropy has 
been studied in [8,9]. Covariant proposal for computing entanglement entropy was formulated 
in [10]. In [11] entanglement entropy in Little String Theory has been studied. Other interesting 
work includes [12-14]. 
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authors interpreted this phase transition as a signature of confinement. In the case of a 
disk in 2+1 dimensions or a ball in 3+1 dimensions there is again a single parameter R, 
and as we will see below there are again two types of solutions in the bulk. The analog 
of a connected solution has a disk topology and dominates at smaller values of R. In this 
solution the circle (sphere) which is the boundary of the disk (ball) shrinks to zero size in 
the bulk, while Kaluza-Klein circle remains finite. The analog of the disconnected solution 
has a cylinder topology and dominates at larger values of R. This is the solution where 
the Kaluza-Klein circle shrinks to zero size. For intermediate values of R the structure has 
a classic "swallowtail" shape typical of first order phase transitions. 

To compute the topological entanglement entropy in 2+1 dimensions, we need to 
consider a disk whose radius is very large compared to the correlation length. In this 
regime the only available solution has the topology of a cylinder. We show that this 
solution approaches the straight cylinder as -R^oo, and the topological entropy associated 
to it vanishes. In hindsight, this is not very surprising, since we do not expect long-range 
topological order in the ground state of 2 + 1 dimensional QCD. Our calculation hence 
can be viewed as a consistency check on both the holographic prescription for computing 
entanglement entropy and topological entanglement entropy. Note that the existence of 
phase transition between the two topologies is crucial for the whole picture to be consistent. 

We then compute the entanglement entropy in the theory on DA branes compactified 
on a circle (this theory has many features of QCD in 3+1 dimensions, as emphasized in 
[16]). The structure of the solutions is similar to the 2+1 dimensional case. There is again 
a first order phase transition between the hypersurfaces of disk and cylinder topology, as R 
is varied. However when computing the radius-independent term we encounter a surprise: 
it contains an explicit dependence on the cutoff. We attribute this to the unconventional 
UV properties of the theory (see discussion Section). 

The rest of the paper is organized as follows. In the next Section we review the 
definitions of entanglement and topological entanglement entropies and the holographic 
proposal of Ryu and Takayanagi. In Section 3 we start with a warm-up example of a 
cylinder in M = A super Yang Mills in 3+1 dimensions and investigate the structure of 
the solutions. We then compactify the theory on a circle to generate a confining theory in 
2 + 1 dimensions. We re-investigate the structure of the solutions, compute the value of 
entanglement entropy as a function of R and extract the value of topological entanglement 
entropy, which vanishes. In Section 4 we repeat this analysis for the case of the DA branes 
compactified on a circle. We find a very similar phase transition structure, while the 
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generalization of topological entanglement entropy turns out to be cutoff-dependent. We 
discuss our results in Section 4. Appendix contains some details regarding the structure 
of the minimal hyper surfaces. 

2. Review of the background material 

2.1. Entanglement entropy and the holographic prescription 

Consider a pure quantum state in a system that can be subdivided into two 
subsystems A and B. If we trace its density matrix over the degrees of freedom of £?, 

PA = Tr B mW > (2-1) 
the density matrix p A will be in general in a mixed state. Its von Neumann entropy 

S A = -Tr A (p A lnp A ) (2.2) 

is called the entanglement entropy and is a measure of the entanglement between A and B 
in the original state |\&). In particular, it is zero if |\I>) is the product of a state in A and 
a state in £?, and one can verify that S a = Sb- 

Of particular interest is the case when the two subsystems A and B correspond to 
two regions in the spacetime of a local quantum field theory and is the ground state 
[17,18]. In this case the entanglement entropy is cutoff dependent and its leading term is 
proportional to the area of the surface that separates the two regions. 

For quantum field theories which have gravitational duals [19], the authors of [5] 
proposed a simple geometric prescription for computing the entanglement entropy of the 
vacuum. For a d + 1 dimensional conformal field theory with an AdSd+2 dual, the idea is 
to find a d-dimensional surface Y which minimizes the action 

and approaches the boundary of the regions A and B at the boundary of the AdSa+2 
manifold. The surface Y is defined at a fixed time and G^ d is the induced string frame 
metric. As shown in [5], for the case of AdS%, this prescription reproduces precisely the 
known expression for the entaglement entropy in 2D CFT, and has been subjected to 
several tests. 

When the boundary theory is not conformal, which is the case we will be interested 
in this work, the above functional generalizes to 

^ = i^0)// Ve " 2 *V / GE (2-4) 

where now Y extends in all remaining directions. 
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2.2. Measures of quantum order and topological entanglement entropy 

Following the discovery of the Fractional Quantum Hall states [20], it became clear 
that there are quantum systems in which different phases or orders (e.g. quantum Hall 
states with different filling fractions) preserve the same symmetries. At zero temperature, 
different orders have different quantum correlations between the microscopic degrees of 
freedom, but cannot be distinguished by any local order parameter. For these quantum 
orders, Landau's theory of symmetry-breaking and local order parameter is inadequate, 
and novel quantum numbers are needed to characterize them. 

When the system has an energy gap the quantum order is called topological order [21], 
and the theory is described in the infrared by a topological quantum field theory. Useful 
characterizations of topological orders are the degeneracy (in the infinite volume limit) of 
the ground state in a Riemann surface as a function of the genus [22,1], the spectrum of 
quasiparticle excitations [23] and the structure of the gapless edge excitations [2]. 

Another quantum number which characterizes topological order is the total quantum 
dimension, defined as follows. The number of linearly-independent states having A quasi- 
particles of type a, for large A, is proportional to , where d a is the quantum dimension 
of the quasiparticle a [24] . The total quantum dimension of the system is 



In general, T> > 1 signals a topological order. In theories with no quasiparticle excitations, 
as in the model we consider in the next section, only the identity sector contributes to (2.5), 
so we expect T> = di = 1. 

As it turns out, the total quantum dimension D is intimately related to the entangle- 
ment properties of the ground state. Consider, in 2 + 1 dimensions, a disk A with a smooth 
boundary of length L in the infinite plane. It was shown in [3,4] that the entanglement 
entropy of the ground state between the disk A and its exterior behaves as 



where a is non-universal and cutoff dependent and the additional terms vanish in the 
limit L — > oo. The quantity 




(2.5) 



S A = aL-log£>H 



(2.6) 



7 = log V 



(2.7) 
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is called topological entanglement entropy. It is a measure of topological order encoded 
in the wave-function, and not in the spectral properties of the Hamiltonian, as the other 
characterizations of topological order mentioned above. For theories in 2 + 1 dimensions 
whose boundary degrees of freedom live in a 2D conformal field theory, the topological 
entanglement entropy can be expressed in terms of the modular matrix of the latter. For 
studies exploiting this idea see [25,26]. 

The definition of 7 can be naturally generalized to dimensions higher than 2 + 1 by 
considering a ball instead of a disk, but these cases have not been thoroughly studied yet. 2 



3. D3 branes on a circle 



Here we consider entanglement entropy for a disk in the strongly coupled four-dimensional 
M = 4 SYM compactified on a circle. At low energies this theory reduces to the three- 
dimensional gauge theory with confinement, mass gap and finite correlation length. The 
metric can be written as 



dU 2 
h(U) 



+ dx^dx^ + h{U)dx\ 



+ L 2 dnl ^ = 0,1,2 (3.1) 



where 



u 2 - — — 

U ° ~ 4R 2 



(3.2) 



and R3 is the radius of the Kaluza-Klein circle, and L 4 = 4tt\ = 47rg s N c sets the curvature 
scale of the AdS space. 3 It is convenient to switch to the variable z = L 2 /U. The metric 
(3.1) takes the following form 



ds' 



dz' 



z 2 h(z) 



+ 



' ix >^ + h{z) ^i + d nl 



where 



h(z) = 1 - 



z = 2R 3 



(3.3) 



(3.4) 



We follow [5], where the prescription to compute entanglement entropy was formulated in 
the holographic setup. Consider a disk in the x\ — xi plane bounded by a circle of radius 
.R, denoted by S\ below. The subscript here refers to the angular coordinate in the x\ — x-i 



2 Solvable systems exhibiting topological order in 3 + 1 dimensions were studied in [27]. 

3 Here and in the rest of the paper we set a! = 1. 
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plane; the radial coordinate is denoted by r. As reviewed in the previous Section, the 
entanglement entropy for this disk can be computed by evaluating (2.3) on the minimal 
8-dimensional surface at t = const which asymptotes to x S^. 3 x S 5 on the boundary 
at z = 0. In this computation we need to set = 87r 6 g 2 ■ The minimal surface is 

(8) 

determined by specifying z(r), with the induced metric G\ nd given by 



a6 ind 



1 + dr 2 + + ^$dx 2 + dQl 



z 2 \ h(z) J z 2 z 2 



(3.5) 



3.1. Conformal limit 



We start with the warm-up example, which corresponds to the conformal theory, where 
the bulk geometry is AdS$ x S 5 . In this limit R^—^oo, h(z) = 1 and the region in the 3+1 
dimensional boundary theory is the cylinder whose length we denote by / (of course, in the 
compactified case, / = 2nRs). Eq. (2.3) can be written as 



s = AN2J 



J dr^l + (z') 2 (3.6) 



15tt 

which gives the following equation of motion for the hypersurface 

rz' \ 3rVl + (z') 2 




Vi + (0 2 J z 
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(3.7) 



One needs to specify the limits of integration in (3.6). For the hypersurface which asymp- 
totes to the cylinder at z = defined by x\ + x\ < R 2 , the lower limit is r = ( we 
will elaborate on this in more detail below). Near the boundary the solution of interest 
behaves like z(r)^0 as r^R and the integral in (3.6) diverges near r = R and needs 
to be regularized. We will introduce the upper limit of integration r a by requiring that 
z(r a ) = a, where a— >0 stands for short distance cutoff. In general, entanglement entropy 
is sensitive to the short-distant modes localized near the border of the region, and hence 
is cutoff-dependent. We will need to understand the behavior of z(r) as z^O (and r-^R). 
The leading term is easy to get from (3.7), 

z ~ 2v A Rv A R 3 r (3.8) 

To go further in the expansion, it is convenient to introduce x = 1 — R/r and make a 
substitution 

z = 2R^f(x) (3.9) 
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0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Fig 1. z(r) for R = 0.5, f = 7 x 10 -7 and R = 1, £ = -6.1 x 10 -6 /&Zue curve]; 
R = 1, £ = —2 x 10 -6 [green curve] and z(r = 0) = 1.29, z'(r = 0) = /red curve]. 



which brings (3.7) into the form 

4x[x(l -x)f"(x)-4x 2 f (x) 3 +(4-5x)f (x)]f(x)+f(x) 2 [2-4x-24x 2 f (x) 2 ] 
+ 3(l-x)x(l+4xf'(x) 2 )-12xf(x) 3 f'(x) - 2/(x) 4 = 



(3.10) 



This equation has a solution in terms of the expansion 

f(x) = 1 + a\x + a2X log(x) + b\x 2 + \>2X 2 log(x) + 63a; 2 log(x) 2 + . . . (3-H) 

where 02 = 1/8 and other coefficients determined by the value of a±. Varying the value 
of a± we can vary the trajectory of the minimal surface in the AdS space. In practice, we 
specify initial conditions at the cutoff z = a = 10 -3 in the form 

2 /? 

z(R)=a; z'(R) = (1+0 (3.12) 

a 

Varying £ effectively varies the trajectory z(r) originating at z(R) = 0. A generic trajectory 
looks like the blue curve in Fig. 1. It reaches a maximal value for z, then turns down and 
goes back to the boundary. The corresponding hypersurface asymptotes to two cylinders 
in the boundary theory. Some particular value of £ corresponds to the surface which 
asymptotes to the cylinder of radius R at the boundary (this is the red curve in Fig.l, it 
can also be obtained by requiring z'ir = 0) = 0). There are, of course, curves with smaller 
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Fig 2. —7 as a function of R — r a for R = 1 

values of £. They look like the blue curve in Fig. 2, with the rescaled value of the radius, 
r{z)— >2r(z), and asymptote to two cylinders on the boundary. 

In the following we will be interested in computing the value of entanglement entropy for 
the solution that asymptotes to the cylinder. It is a simple exercise to extract the UV- 
divergent term (see also [7]). For this purpose the parametrization r(z) is more convenient. 
The value of expression (2.3) in this parametrization is 

T2l 



where dot denotes derivative with respect to z and a is the UV cutoff. Inverting (3.9) and 
(3.11) we obtain 



r = R 



z 2 



l-^ + 0(zHogz) 



(3.14) 



this expression can be substituted into (3.13) to obtain 



S = t£{^ + M 1o Zr)-^ (3 - 15) 

where 7 is finite in the limit when short distance cutoff is taken to zero, a^0. Numerical 
evaluation of the integral (3.13) and subtraction of divergent terms in (3.15) produces a 
finite value of 7 = 0.305. In Fig. 2 we plot the (normalized) value of the integral in (3.6) 
and observe that it is indeed finite in the limit a^0, r a ^R. It is interesting that the 
limiting value of 7 is positive and non-zero, similarly to the topological entropy in (2.6). 
Of course, in the conformal theory the correlation length is infinite, so it is not completely 
clear what the physical interpretation of nonvanishing 7 is. 
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3.2. Finite compactification radius 

Consider now the case with finite Kaluza-Klein radius R3, and hence finite correlation 
length. For the rest of this Section it will be convenient to pass to rescaled variables, 

r^ — ; z-> — ; — (3.16) 

Zo Zq Zo 

It is useful to remember that z = 2R 3 defines the scale of the correlation length in the 
theory. The rescaled coordinate z is now bounded from above by z = 1, where the Kaluza- 
Klein circle shrinks to zero size. Eq. (2.3) becomes 



/ Z ° dr^v^T + W (3-17) 
Jo z 



15 

The equation of motion which follows from (3.17), 

= r(z 4 -3[l + (z') 2 }) 
' z^^/1- z 4 + (z') 2 




(3.18) 



in general has to be solved numerically. A detailed study of possible solutions to (3.18) is 
relegated to the Appendix. Here we briefly summarize the results. There are two types of 
solutions which asymptote to the circle of radius R in 2+1 dimensional boundary theory 
(they are shown in Fig. 3). One type involves a connected surface which approaches 
the circle of radius R near the boundary z = and has the disk topology in the (r, 0) 
coordinates. It corresponds to the blue curve in Fig. 3. This solution can be found be 
starting at r = with boundary conditions = z(r = 0) < 1, z'(0) = and integrating 
(3.18) outwards in r. In Fig. 4 the value of R is plotted as a function of z*. Similar to the 
slab geometry, studied in [15], there appears to be a maximal value of R, beyond which 
this solution does not exist. Moreover, in the first approximation there are two branches 
for sufficiently large values of R. The new feature in comparison with [15] is the existense 
of the lower bound on R for the branch which consists of the solutions passing close to the 
end of the space at z = 1. We will call the solution with z*— >1 critical, since, as we will 
see below, it joins with the solution which has the topology of the cylinder, (see below). 
We denote the value of R for such critical solution by R c . 

Another interesting solution, which was absent in the conformal case involves the 
surface which starts at z = 1, r = tq and goes all the way to the boundary. This surface has 
the topology of the cylinder in the (r, 4>) coordinates and is an analog of the disconnected 
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solution in [15]. It corresponds to the red curve in Fig. 3. To exhibit this solution (and 
for the purpose of computing the entanglement entropy below) it is convenient to use 
parametrization r(z), like in (3.13). In the non-conformal case the action is given by 



f ■ v 

s = 



and the equation of motion 



d (1 - z^)rf _ y/j + (1 - z 4 )r 2 
d~z z 3 v^TTT - z 4 )r 2 ~ ~ 3 

Near z = 1 (3.20) has a solution 



(3.20) 



r = r + V^ + --- ( 3 - 21 ) 

where ro is the value of r at z = 1 and the dots stand for the terms subleading in 1 — z. 
Starting with these boundary conditions and integrating (3.20) or (3.18) numerically we 
obtain the solution with the topology of the cylinder. Note that at large r the solution of 
(3.21) can be found in the form 

r = r + ^ + o(l/r ) (3.22) 

We can also verify (3.22) numerically. In Fig. 5 we show the behavior of R as a function 
of r . At large r we recover (3.22), but at small r there is some interesting structure. 
As r — »0, the value of R approaches R c and the solution approaches the critical solution 
(with the disk topology) discussed earlier. As ro increases, R goes down, before climbing 
up back again. 

In addition to the two solutions described below, and just as in the conformal case 
studied above, there are also solutions which approach the two concentric circles on the 
boundary. Such solutions will define an annulus in the boundary 2+1 dimensional theory. 
We will not discuss them further; more details are provided in the appendix. 
To compute the entanglement entropy we follow the steps of the previous subsection, 
where the conformal case was discussed. In particular, the divergent terms are again given 
by (3.15). We define the normalized finite part of the entropy S via 

15 V « 2 4# B RJ 15 V ; 

10 



0.5 1.0 1.5 2.0 



Fig 3. z{r). Disk [blue] and cylinder [red] topology. 
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Fig 4. R as a function of z* (solutions of disk topology) 

The value of S as a function of R [normalized according to (3.16)] is plotted in Fig. 6. 
The blue part of the curve corresponds to the solution which has the disk topology [blue 
curve in Fig. 3]. It is (modulo some fine details) double- valued, since R is a double- valued 
function of be inferred from Fig. 4. For a given value of R the branch with 

smaller z* dominates the computation of entanglement entropy. The solution with z*— >1 
smoothly connects to the solution which has the topology of the cylinder (red curve in 
Figs. 3,6). The value of S for this solution is again double- valued, according to Fig. 5. 
As the value of R is increased, the cylinder solutions begins to dominate the entanglement 
entropy (this happens at the point where the red line crosses the blue in Fig. 6) . For larger 
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Fig 6. S as a function of R; 

values of R only the solution with cylinder topology exists. According to (3.22), for large 
R, 

r = R + 0(R~ 1 ) (3.24) 

and hence 

S = -- + 0(R~ 1 ) (3.25) 
2 

which implies vanishing topological entropy 7. 
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4. D4 branes on a circle 



In this section we will consider the entanglement entropy of a three-dimensional ball 
in a confining theory in d = 3 + 1. The theory is obtained from compactifying a 4 + 1 
theory with radius R4 [16]. At strong coupling, it has a gravitational description with 
metric, 



-20 



3/2 



3/2 



R\ 3 dU 2 



UJ f(U) 



+ dx^dx^ + f{U)dx\ 



+ R 3 / 2 U 1/2 dfli f/ = 0,l,2,3 



(4.1) 



where x± ~ X4 + lnR^ and 

f(U) = 1 - 
The constant R is given by 



where 



Uo 
U 



R s = ttA, 



A = (27r)- 2 AT c ^ M = g s N c l s 



U 



4n_\_ 
"9" if 



(4.2) 



(4.3) 



is the five-dimensional 't Hooft coupling. The supergravity description is valid as long 

as < 1, 

U<^R. (4.4) 

The boundary is at U — > 00, and it is convenient to change the (U, x^) coordinates 
into (z, x^) defined by 



U 



Uo 



x h 



Uo 



3\ 1/2 



(4.5) 



x h 



We have now z G [0, 1] and the boundary is at z = 0. The metric and dilaton are now 

dz 2 



ds 2 = (R 3 U ) 1/2 I z~ 3 / 2 



+ dx^dx^ + (1 - z 3 )dx\ 



z(l-z 3 ) 



+ z~ 1 / 2 dn 2 4 



-20 



R 
U~o 



3/2 



(4.6) 



,3/2 
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We consider dividing the 3 spatial directions into a three-dimensional ball and its 
exterior, such that the surface between the two domains is a two-sphere. Therefore, we 
are interested in an 8-dimensional surface in the bulk, which wraps the inner 4-sphere, 
is constant along the periodic direction x 4 and ends in a two-sphere of radius p in the 
boundary z = 0. 

Let us call r the radial coordinate in the (x 1 ,x 2 ,x 3 ) directions. The 8-dim surface 
is specified by the curve z(r), or alternatively, by r(z). Choosing z as the independent 
variable, the induced metric is 



ds 2 = (R 3 U ) 1/2 {z~ 3 / 2 



1 



and the area functional (2.4) is 



+ f z {z) ) dz 2 + r 2 (z)dn z 2 + (1 - z 6 )dx'i 



3\ J~2 



+ z~ 1/2 dQ 2 4 
(4.7) 



S 



1 



/ 



(10) I a ae V ^ind 



n s R 6 U R4 



2G 



(10) 
N 



dz 



r 2 (z) 



y/l + z(l - z 3 )r> 2 (z) 



(4.8) 



where we have used Vol(S 2 ) = ^ and Vol(S 4 ) = ^p. The equation of motion from this 



action is 



4-4z(-l + z 3 ) r'[z l2 
+ r[z) (2 (2 + z 3 ) r'[z] +z(5-7z 3 + 2z 6 ) r'[zf + 2z (-1 + z 3 ) r"[z]) = 



(4.9) 



Alternatively, choosing r as the independent variable, the induced metric and the area 
functional are 



ds 2 = (R 3 U ) 1/2 \z- 3 / 2 

n 3 R 6 U R4 



1 + f (r \, ) dr 2 + r 2 dQ 
z[l — z 6 ) J 



2G 



(10) 
JV 



l -f (1 - z 3 )dx\ + z- 1/2 dill 

(4.10) 

f dr^—-^z(r)(l - z 3 (r)) + z' 2 (r) (4.11) 
J z°{r) 



and the equations of motion are 



2rz [r] 7 + 2rz [r] 3 z' [r] 2 + Az' [r] 2 (r + z' [r] ) 

+ z[r] {Az'[r} + r (5 + 2z"[r))) - z[r} 4 (4z'[r] + r (7 + 2z"[r])) = 



(4.12) 
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4-1. The two topologies 

The shape of the surfaces we are looking for is determined by the functions r(z) or z(r). 
As in the previous Section, starting with z = at r = p, there are two possible shapes for 
z(r): 

1) the curve starts from z = at r = p and reaches z = 1 at some value r = tq. We 
call these curves of 'cylinder' topology. 

2) the curve starts from z = at r = p and ends at r = with z = zq < 1. To avoid 
a conical singularity, we should demand ^'(0) = 0. We call these curves of 'disk' topology. 

In practice we can obtain the curves numerically. For the disk topology, the initial 
conditions of the numerical problem are z(0) = z$, z'(0) = 0. 4 

For the cylinder topology, the numerical problem can be set to find either z(r) or r(z), 
and in both cases we can set the boundary conditions either at z = or at z = 1. Consider 
the latter case first. With the boundary condition z(r ) = 1, the leading term in the 
solution of (4.12) is 

3 

z{r) ~ 1 - -r (r - r ) (4.13) 
One can check that inverting this function, 

2 

r{z) ~r + — (1-z), (4-14) 
Sr 

gives the leading term of the solution of (4.9) with b.c. r(l) = r . From the expres- 
sions (4.13) or (4.14) one can read out the boundary values of z'(r ) or r'(l) needed for 
the numerical solution. 

Alternatively, we can consider the behavior of the solutions near z = 0. The leading 
term in the solution of eq.(4.12) near the boundary (z ~ 0), with boundary condition 
z(p) = 0, is 5 

z(r) ~ p(p — r) + ■ ■ ■ (4-15) 



4 One can also try to find numerically the 'disk' solutions by fixing the boundary conditions 
at z = (instead of at r = 0), but the shape of the curve is very sensitive to the initial conditions 
and a slightly bigger or smaller value for z'(p) does not give a solution ending at r = with 
z'(0) = 0. 

5 There is a second solution to (4.12) that behaves as z(r) ~ — \{p — r) 2 , but it is not physical 
since z € [0, 1]. 
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The inverse to the above function has the expansion 

which are the leading terms of the solution of (4.9) with b.c. r(0) = p. One can verify 
that the series (4.16) contains terms of the type z n In 771 z, starting with z^lnz. From the 
expressions (4.15) or (4.16) one can read the boundary values of z'(p) or r'(0) which can 
be used in the numerical solutions. 

In Figure 7 we show the profile of z(r) for solutions of both topologies. As zq ap- 
proaches 1, there is a smooth transition between the two topologies, since the 'disk' solution 
with zq = 1 coincides with the 'cylinder' solution with ro = 0. 




The value of r(0) = p has an interesting dependance on the value of zq and tq for the 
'cylinder' and 'disk' solutions, respectively, as we show in Figs. 8 and 9. 

We see that for the 'disk' solutions there is a maximal value p as zq grows, and for 
bigger zq the value of p shrinks. A more detailed numerical study, which is not evident in 
Fig. 8, shows that the value of p actually oscillates with decreasing amplitude as Zq — > 1. 
For the 'cylinder' solutions, p grows monotonically with ro for large tq. As tq diminishes, 
there is a minimum value for p, such that for smaller values of ro, p starts to grow again. 
Also here a close look shows that the value of p oscillates with decreasing amplitude 
as ro — > 0. 
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Fig 8 p as function of zq for the 'disk' solutions. 




Fig 9 p as a function of r$ for the 'cylinder' solutions. 

4-2. The phase transition 

The entanglement entropy is UV- divergent and we will be interested in subtracting 
the divergencies and obtaining a finite quantity. The divergent terms in the entanglement 
entropy can be obtained by inserting the leading behavior (4.16) into (4.8). This gives 



S = 



7r 3 R 6 U Q R 4 



2G 

r3D6 



(10) 
N 



dz 



r 2 (z) 



z- 



'l + z{l-z*)r' 2 (z) 



7T R U0R4 



2G 



(10) 
N 



2a 2 



3 In a 
2a + V 



+ S 



finite 



(a) 



(4.17) 
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where Sfi n i te {o) is finite as a — > 0. 6 Unlike the slab geometry, the divergent term does 
depend on the size p. In Fig. 10 we show the value Semite (0) as a function of p. Note that 
this is a multivalued function, since a given p can correspond to two different values of z 
or ro, as we saw above. 




u 1.475 1.480 1.485 1.490 1.495 

Fig 10 The finite entanglement entropy 5'/mite(0) as a function of p. The red line 
corresponds to the 'disk' topology and the blue line to the 'cylinder' topology. 



4-3. Topological entropy 

To compute the topological entropy only the cylinder solutions are relevant, since we 
are interested in the the limit p — > oo. For the same reason, it is convenient to expand the 
solution in powers of p as 

r(z)=p + f(z) + 9 -^ + ^ + --- (4.18) 
P P 

where /(0) = g(0) = h(0) = . . . = 0. Inserting this expansion in the equation of motion 
(4.9) for r(z), gives for the leading term (~ p), 

4f'[z] + 2z 3 f'[z] + 5zf'[z] 3 - 7z 4 f'[z] 3 + 2z 7 f'[z] 3 - 2zf"[z] + 2z 4 f"[z] = (4.19) 



6 It is interesting that from the r 2 and r' 2 dependence in (4.17) there are contributions from 
the z 2 term in (4.16) to the In a term in (4.17), but one can check that they cancel and the 
divergence would be the same had we kept only the linear term in (4.16). 
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This equation is solved by f'(z) = 0, and from /(0) = we get that f{z) = 0. With this 
result, the next term in the equation of motion (4.9) for r(z) (~ p°) is 



4 + 4g'[z] + 2z 3 g'[z] - 2zg"[z] + 2z 4 g"[z] = 
which has the exact solution 



1 + cz^ 



(4.20) 



(4.21) 



1-z 3 ' 

and c is an integration constant. The 'cylinder' solutions are required to be regular at z = 1, 
and this fixes c = — 1, so we have 

1 + z 



9\z) = -- 



+ Z + Z" 



(4.22) 



Integrating this function gives 

arctan 
9{z) = — 



1+22 

V5 



V3 



1 log [1 + z + z 2 ] + 

2 L J 6\/3 



, z 4 z 7 . ,z 3 z % z 9 . 



(4.23) 



where the integration constant is fixed by g(0) = 0. 

As it turns out, the 'topological' term in the entanglement entropy depends only 
on g(z) and therefore can be computed exactly. Indeed, inserting the expansion (4.18) 
into (4.8) gives, using /(0) = 0, 



tt 3 R 6 U R 4 



S = 



2G 



(10) 
N 



dz 



p 2 4g[z]+z(l-z 3 )g'[z]'- 



z 6 



= Sp ,- 1 + 0(p- 1 ) 



where 



n 3 R 6 U Q R 4 2 / 1 



S p 2 = 



7 = 



AG 



(10) 
N 



-^-1 = 



0/ 



2^ 3 



n 3 R 3 R 4 r 2 



+ 0(p~ 1 ) 



(4.24) 



4G 



(10) 
N 



tt 3 R 6 U R 4 fir 



8G 



(10) 
N 



+ 2 + log3- - 
3 a 



(4.25) 



7T 3 R 6 U R 4 (it n , „ 6C/oo 



8G 



(10) 



where C/qq = C/o/a, and the dimensionful radius of the ball is = p 2 R 3 /U . In order to 
express the entropy using gauge theory quantities we should use 

8tt 6 A 2 



gT = 8A2 = 



N 2 



(4.26) 
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and the relation between Uqo and the field theory cutoff A in our background [28] , 

Uo, = 9ymN c A 2 = (2tt) 2 AA 2 (4.27) 



This gives finally 



S p2 = XNtr'IU 



2 2 . / TT'-'A' 



(4.28) 

So we find a non-zero and cutoff dependent topological entanglement entropy. But we 
should remember that the supergravity description can only be trusted in the limit (4.4), 
which for the cutoff A means, 

A<< (A^' < 4 ' 29 > 



5. Discussion 

In this paper we computed the entanglement entropy for a disk in 2+1 dimensions and 
a ball in 3+1 dimensions in the confining theories with holographic duals. We observed 
a phase structure similar to the one found in [15], with two types of solutions. At small 
values of R the solution of the disk topology dominates the computation of entanglement 
entropy, while at large values of R the cylinder type solution dominates. In 2+1 dimen- 
sions this structure ensures vanishing topological entropy, which is an expected result for a 
QCD-type theory. This can also be viewed as a consistency check of the holographic pre- 
scription. (Of course, it would be great to make progress in understanding the prescription 
of Takayanagi and Ryu from the first principles.) It is interesting that in the conformal 
case the analog of the topological entropy is nonzero. It is not completely clear what the 
physical interpretation of this fact is. 

In 3+1 dimensions one can define a quantity which is analogous to topological entropy. 
It turns out to be UV divergent. A possible explanation of this behavior is that, as 
mentioned, the supergravity description of the field theory cannot be trusted beyond the 
bound (4.29). On the other hand, the topological entropy has not been studied much in 
3+1 dimensions, and such studies might shed light on our result. 

Both in the 2+1 and 3+1 cases the solution with the disk topology joins the solution 
with the cylinder topology at some critical value of the disk or ball radius. The behavior of 
the entanglement entropy near this point seems reminiscent of self-similar type of behavior 
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of free energy of probe branes, described in [29] . It would be interesting to investigate this 
in more detail. 

For the slab geometry, signatures of the phase transition predicted in [15] have recently 
been found in numerical studies in the lattice [30] (see also [31]). It would be interesting 
to seek for similar signatures of the phase transition predicted by the disk and spherical 
geometries. But it should be noted that, unlike the case of the slab geometry [15], in our 
case the (regularized) entropy scales as AT? at both sides of the critical point. This suggests 
that the existence of the phase transition might be unrelated to the confining properties of 
the theory. Studying entanglement entropy for other geometries, different from disks and 
slabs, might also shed light on the relation of entanglement entropy to the properties of 
the theory. 

Finally, it would be interesting to find supergravity backgrounds which would cor- 
respond to theories with non-zero topological entanglement entropy. One instance could 
be the system of D2-D6 branes which models the Quantum Hall effect [32] , but to apply 
the methods of this paper one would need the full supergravity solution, and not only the 
probe approximation used in [32]. 

Acknowledgements: We thank S. Giombi, C. Herzog, I. Klich, M. Levin, M. Kulaxizi 
and D. Kutasov for discussions. Ari Pakman thanks Harvard University for hospitality. 
The work of Ari Pakman is supported by the Simons Foundation. Andrei Parnachev thanks 
KITP, Santa Barbara, University of Amsterdam, Ecole Polytechnique and University of 
Washington for hospitality during the course of this work. 

Appendix A. Solutions in detail 

To investigate the structure of the solutions of (3.18) and (3.20) we repeat the analysis 
of Section 3.1. The UV behavior of the solution of (3.18) still satisfies (3.9) and (3.11), 
although eq. (3.10) is modified. Typical behavior of the solutions for sufficiently small 
R are shown in Fig. 11. The picture is completely similar to the one in Fig. 1, which 
is not surprising, since R is substantially smaller than the correlation length. In Fig. 
12 r(z) is computed. In particular, there is a maximal value of z which the curves can 
attain. The red and blue curves in Fig. 12 approximate to the annuli with internal radius 
Ri = R = 0.468 and external radius R e > R. The green curve in Fig. 12. (which is the 
counterpart of the green curve in Fig. 11) gets very close to the connected solution but 
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eventually approximates to the annuli with R e = R = 0.468, Ri <C R. In Fig. 13 we plot 
the structure of the solutions for the critical value of R = R c = 0.7034. The blue curve 
is essentially the critical solution, while the red curve comes very close to z = zq before 
turning up and ending back on the boundary again. It is instructive to see how this picture 
changes for R > R c in Fig. 14. We now start with R = 1 > R c . The red line gives rise to 
the annulus with R4 < R e = 1, while the blue line corresponds to R e > R4 = 1 and comes 
very close to z = zq. The green line separates the two domains; it is natural to think that 
it ends at z = z$ and gives rise to a surface with annulus topology in the bulk (not to be 
confused with the topology of the boundary which is a circle). This is the analog of the 
disconnected solution in [15]. This solution can be exhibited in a more direct way. One 
can analyze the equation of motion near z = l,r = tq > and integrate the solution from 
there. This is done in Section 3.2. 




0.2 0.3 0.4 

Fig. 11 z(r) in confining background for R = 0.468; £ = —5 x 10 -4 [blue curve]; 
£ = —8.05 x 10 -6 [green curve]; z'[r = 0) =0 [red curve] 
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